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Abstract 

In this paper we present a robust mixed integer optimization model to utilize re¬ 
generative braking energy produced by trains in a railway network. An electric train 
produces regenerative energy during braking, which is often lost in present technology. 

To utilize this energy we calculate a timetable which maximizes the total overlapping 
time between the braking and accelerating phases of suitable train pairs, so that the 
regenerative energy of a braking train can be transferred to a suitable accelerating 
one. We apply our optimization model to different instances of a railway network for 
a time horizon spanning six hours. For each instance, our model finds an optimal 
timetable very quickly (largest running time being 86.64 seconds). Compared to the 
existing timetables, we observe signihcant increase in utilization of regenerative energy 
for every instance. 

Index Terms — Railway networks, energy efficiency, regenerative braking, train schedul¬ 
ing, mathematical optimization. 


1 Introduction 

1.1 Background and motivation 

Application of mathematical optimization techniques to efficient energy management of elec¬ 
tric vehicles has gained lot of interest in the past few years[l, 2, 3, 4]. Among all of the public 
transport modes, trains are often preferred by the passengers for providing higher capacity 
and safety. In most railway networks, trains use electricity as their primary source of energy 
and are equipped with regenerative braking technology. When a train makes a trip from an 
origin platform to a destination platform, its optimal speed prohle consists of four phases: 
1) accelerating, 2) speed holding, 3) coasting and 4) braking [5], as shown in Figure 1. in 
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a qualitative manner. Most of the energy required by the train is consumed during the 
accelerating phase. During the speed holding phase the energy consumption is negligible 
compared to accelerating phase, and during the coasting phase there is no need for energy. 
When the train brakes, it produces energy by using a regenerative braking mechanism. This 
energy is called regenerative braking energy. Naturally, proper utilization of regenerative 
braking energy of trains can lead to signihcant energy savings. 

However, transfer of the regenerative braking energy back to the electrical grid requires 
specialized technology such as reversible electrical substations [6, page 30]. On the other 
hand, storing it using current technology is not viable, as the storage options e.g., via super¬ 
capacitors, flywheels [7] are very expensive and have high discharge rates [8, page 66], [9, 
page 92]. A better strategy to utilize the regenerative braking energy of a train, that can 
be implemented with current technology, is to synchronize its braking phase with the ac¬ 
celerating phase of another nearby train operating under the same electrical substation. A 
positive overlapping time that arises from such a synchronization process enables transfer 
of the regenerative braking energy of the hrst train to the second one via the overhead 
contact line or a third rail [10], and can save the electrical energy that would be lost other¬ 
wise. Our objective is to design a railway timetable that contains the arrival and departure 
time of every train to and from all the platforms it visits, such that the duration of the 
synchronization processes between suitable train pairs (SPSTPs) is maximized subject to 
the different constraints present. 

1.2 Related work 

Over the past three decades, the general timetabling problem in a railway network has been 
studied extensively with objectives such as total cost, total passenger delay, quality of the 
service etc. [11]. However, a few among them attempt to determine a timetable to utilize 
regenerative braking energy of trains; all of them have appeared in the last three years. 

The work by Peha-Alcaraz et al. [12] formulates a Mixed Integer Programming (MIP) 
model. In this model the objective is to maximize the total duration of all possible synchro¬ 
nization processes between all train pairs. The model is then successfully applied to line 3 of 
the Madrid underground system in Spain. However, for many railway networks, inclusion of 
all the train pairs in the objective may not be realistic, because energy transfer between trains 
that are not close to each other suffers from drastic transmission loss. Besides, inclusion of 
all such train pairs can make the optimization problem computationally intractable even for 
a medium sized networks, as each such train pair in the model corresponds to one binary 
variable. Because of the absence of connection and turn-around constraints, the model is 
only applicable to a single train-line. 

A more tractable MIP model is presented in [13]. The model considers only the suitable 
train pairs for transfer of regenerative energy and can be applied to railway network of 
arbitrary topology. The optimization model is then applied to the Dockland Light Railway. 
The results show a signihcant increase in the total duration of the synchronization process. 
Although such increase, in principle, may increase the total savings in regenerative energy, 
the actual savings in terms of energy is not calculated. Another drawback of this model is 


2 



Speed 



Figure 1: The speed profile of a train (top), and corresponding energy consumed and produced 
in accelerating phase and braking phase respectively (bottom). 
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that it lacks robustness, as it assumes that even if the trip time changes, the duration of the 
associated braking and accelerating stay the same, which is not the case in real life. 

Other relevant works implement meta-heuristics such as simulated annealing [14] and 
genetic algorithm [14] and [15]. These models cannot provide any guarantee of optimality 
for the generated timetable, because a meta-heuristic approach depends on manual tuning 
of the parameters. 

An analytical study of a simplihed and periodic railway schedule appears in [16] with 
focus on storage of regenerative energy. The model does not incorporate periodic event 
scheduling constraints [17] used to formulate periodicity in a timetable, so it is not possible 
to apply the model to some cases, e.g., where trains associated with the second period enter 
the network while the trains associated with the hrst period are still running. 

1.3 Contributions 

Our contributions in this paper can be summarized as follows: 

• Using mathematical programming techniques, we propose a robust MIP model to uti¬ 
lize the regenerative braking energy of trains in a railway network. In comparison 
with the existing works, our model can be applied to a railway network of arbitrary 
topology. 

• For most of the existing railway networks, the railway management has a feasible 
timetable, which we exploit to devise an optimization strategy. This strategy produces 
a smaller search space and is computationally more tractable in comparison with the 
related works. We prove that all possible cases arising from an SPSTP are modelled 
accurately by our model. 

• We apply our optimization model to different instances of a real-life railway network 
and hnd that our model gives optimal timetables very quickly. We observe signihcant 
increase in the utilization of regenerative braking energy. 

1.4 Organization 

This paper is organized as follows. In Section 2 we describe the notation and notions to be 
used in this paper. In Section 3 we dehne an SPSTP mathematically, propose our optimiza¬ 
tion strategy and formulate a mixed integer linear objective function. Then in Section 4 we 
model and justify various constraints present in the railway network. We propose the full 
optimization model in Section 5. In Section 6 we apply our model to realistic instances of 
different size. Section 7 presents the conclusion. 


2 Notation and notions 

All the sets described in this paper are strictly ordered and hnite unless otherwise specihed. 
The cardinality and the ith element of such a set S is denoted by liSI and S{i) respectively. 
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Consider a railway network where the set of all stations is denoted by S. The set of all 
platforms in the railway network is indicated by Af. A directed arc between two distinct 
and non-opposite platforms is called a track. The set of all tracks is represented by A. The 
directed graph of the railway network is expressed by ^ = (A/*, A). A train-line or line is 
a directed path, where the set of nodes represents non-opposite platforms and the set of 
arcs represents non-opposite tracks. A crossing-over is a special type of directed arc that 
connects two train-lines. If a train arrives at the terminal platform of a train-line, turns 
around by traversing the crossing-over, and starts travelling through another train-line, then 
the same physical train is treated and labelled functionally as two different trains by the 
railway management [18, page 41]. The set of all trains to be considered in our problem is 
denoted by T. The sets of all platforms and all tracks visited by a train t in chronological 
order are denoted by A/”* C A/ and A^ ^ A respectively. The train-path of a train is the 
directed path that contains all platforms and tracks visited by it in chronological order. 

The decision variables to be determined are the arrival and departure times of trains to 
and from the associated platforms respectively. Let a- be the arrival time of train t E T 
at platform i G AT^ and d* be the departure time of train t from platform j G Af^. Our 
objective is to maximize the duration of total overlapping time of all SPSTPs subject to 
the constraints, which would make it possible to save signihcant amount of electrical energy 
produced by the braking trains by transferring it to the accelerating trains. 


3 Modelling the objective 

At hrst, we need to characterize the train pairs and the associated platform pairs necessary 
to describe the SPSTPs. The platform pairs to consider are those opposite to each other 
and powered by the same electrical substations, because the transmission loss in transferring 
electrical energy between them is negligible. The set that contains all such platform pairs is 
denoted by fl. Consider any such platform pair {i,j) E fl, and let 7i C T be the set of all 
the trains which arrive at, dwell and then depart from platform i. Suppose, t Ell. Now, we 
are interested to hnd another train t on platform j, i.e., t E Tj, which along with t would 
form a suitable pair for the transfer of regenerative braking energy. To achieve this, we 
start with an initial feasible timetable for the railway, which represents the desired service 
to be delivered. Too much deviation from it is not desired. For every train t, this timetable 
provides a feasible arrival time a* and a feasible departure time d\ to and from every platform 
i E Af* respectively. Intuitively, among all the trains that go through platform j, the one 
which is temporally closest to t in the initial timetable would be the best candidate to form 
a pair with t. The temporal proximity can be of two types with respect to t, which results 
in the following dehnitions. 

Definition 1. Consider any {i,j) ^ For every train t E the train t E Tj is ealled the 
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temporally closest train to the right of t if 


t = argmin 

of-\-dP 

t'e {xeT^ :0< - ^4^ <r} 


K + d\ d^+d] 


( 1 ) 


where r is an empirieal parameter determined by the timetable designer and is mueh smaller 
than the time horizon of the entire timetable. 

Definition 2. Consider any {i,j) G For every train t ^ the train t ^ Tj is ealled the 

temporally closest train to the left of t if 


t = argmin 

at+d* a?+J¥ 

t'e{a;e7}:0< * * - ^ ^ <?-} 


of+ 4 + 


( 2 ) 


Definition 3. Consider any {i,j) ^ For every train t ^ the train t ^ Tj is ealled the 

temporally closest train to t if 


t = argmin 
t'e{ t,t} 


+ d\ CLj + dj 


( 3 ) 


If both t and t are temporally equidistant from t, we piek one of them arbitrarily. 

Any SPSTP can be described by specifying the corresponding i, j, t and t by using the 
dehnitions above. We construct a set of all the SPSTPs, which we denote by £. Each element 
of this set is a tuple of the form (bj, t, t). Because t is unique for any t in each element of 

£, we can partition £ into two sets denoted by £ and £, which contain elements of the form 
{i,j,t^ t) aiid {i,jd^ t) respectively. For every {i,j,t^ t ) E £ our strategy is to synchronize 
the accelerating phase of t with the braking phase of t. On the other hand, for every t ^ £, 
then it would be convenient to synchronize the accelerating phase of t with the braking 


phase of t. For every {i,j,t^ t ) ^ the corresponding overlapping time is denoted by a. 
and for every {i,jd^ t) G £1, the corresponding overlapping time is denoted by a* 


ij 1 

Our 


objective is to maximize the sum of overlapping times over all the elements of £ and £, i.e., 


maximize 


a, 


tt 

ij 


E + E _ 

subject to the constraints present in the system. 


( 4 ) 


We model aC for all {i,jd, t) E £ and aC for all {i,jd, t) E £ in terms of the arrival 
and departure times of trains. Consider the case, when {i,jd^ t ) E £. We need to ensure 


6 



that after we apply the optimization strategy, t still stays the temporally closest train to the 
right of t. Otherwise, the only way to achieve a positive overlapping time is to synchronize the 
braking phase of t with the accelerating phase of t , which might result in a large deviation 
of event times compared to the original timetable, especially when there is no or very little 
overlapping to begin with. We write this constraint as follows: 

fa/ + dj - a\- d{) 

\ > 0 . ( 5 ) 

fa/ + d/ - a* - d* + ej 

Here e is a very small positive number to prevent division by zero. 

Let us denote the start of the braking phase of train t before arriving at platform i by 
a*~ and the end of its accelerating phase after departing from the same platform by d-'*'. For 
all t G T and for all i G M*, the duration of the associated braking phase is /?/ = a- — a*“ 
and the duration of the associated accelerating phase is aj = dl'^ — d\. Depending on the 
trip time of the associated trip in consideration, both the durations are within some time 
bounds, i.e; a\ G [a*, a*] and /?/ G [/5*,/3j. Though we do not know the optimal trip time of 
the trains in advance, these lower and upper bounds can be calculated by existing software 
[19, page 3]. For the same reason, the start of the braking phase and end of the accelerating 
phase are within time bounds described by a[“ G [o* —a* —/5*] and dd^ G [d* + a-, d* + a*]. 
All the time bounds are on the order of seconds, as the trip time variation are on the order 
of seconds, so it is reasonable to pursue a robust formulation. To model the overlapping time 

for all t ) E £, we propose the following lemma. 

Lemma 1. For all (i,j, t, t) E £, the overlapping time a// between the aeeelerating phase of 

t on platform i and the braking phase of t on platform j, where {i,j) E D, can be modelled 
by the following robust formulation 


«/- 

- d\ + e < of + 

+ M(1 — A-/ ), 

(6) 

4 - 

-a/ + e < M(1 - 

A“), 

(7) 

< 

IV 

o 


(8) 

^ij 

< , 


(9) 

^ij 

<VA[/, 


(10) 

^ij 

<d\ — aj + a ■ + 

+ 

1 

(11) 

^ij 

^ a/ — d[ + M (1 

-A^). 

(12) 


where M is a large positive number, e is a small positive number smaller than time granularity 
considered and A*/ is a binary variable which is one if and only if is positive. 
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Figure 2: All possible overlapping times between the accelerating phase of train t and braking 
phase of its temporally closest train 8 


























































Proof. We use the hypograph approach to model the overlapping time ajj in terms of the 
associated event times [20, page 75, 134]. Consider any G [a^ — I3\,a\ — /5*] and any 
d!f^ G [d\ + gf, d\ + af\ for some trip time. From Interval algebra [21], we know that there can 
be thirteen different kinds of overlapping possible between the accelerating phase of train t 
and the braking phase of train t as shown in Figure 2. However, there is structure among 
these thirteen relationships when we are concerned with the associated overlapping time aC . 

Among them, when > d*'*' or aj < d\ there will be no overlapping, i.e., aC = 0. We 
will model this scenario using the binary indicator variable A**'. If any of the conditions 

> dff^ or a/ < df occurs, then Af/ is zero, i.e., 

(laf-> rf‘+) V (aj < d‘)) =. (A*/ = 0) (13) 

If is zero, then the overlapping time will be zero, i.e., 

(A“ = 0) ^ (alj = 0) 

> 0) ^ (A*/ = 1), [as a*/ 0]. (14) 

The Equations (13)-(14) are not in mathematical programming format. Using integer 
programming modelling rules [22, pages 166, 172-174, 183-184], we can model them as follows: 


4+ + e< M(1 - A“ ), 

(15) 

ttj + e < M{1 — Xlj ), 

(16) 

< >0, 

(17) 

<MA[/. 

(18) 


When (Oj < and (aJ > d*), we can see from Figure 2 that the overlapping time is 
min{d(''', a| } — max{d*, ~} > 0. Now we want to model the following logical condition 

(A-/ = 1) ^ = min{d[+, a/ } - max{d*, a/ “}) 

using mathematical programming. Using integer programming modelling rules we can model 
this situation as follows: 

t ^ ^t\tt 
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(19) 

( 20 ) 
( 21 ) 



(22) 


< So/ -d\ + Mil-\\‘). 

Combining Equations (15)-(18) and Equations (19)-(22), and using the fact that M ^ 

max{a*, P- }, for any G [a\ — — /?*] and any G [d\ + a*, d- + a*] we arrive at the 

following set of equations: 

a;--d^+ + e<M(l-A“), 
dl — Qj + e < M(1 — \lj ), 

< > 0 , 

< < , 

< <dr-a/- + M(l-A*/), 

<a/ -d^ + M(l-A“). 

Using a*“ G [a* - ^*, 0 * - and d*+ G [d* + a-, d* + a*], the equations above can be 
transformed into the following robust formulation: 

a/ - d^ + e < + M(1 - A“ ), 

d\ — a j + e < M(1 — A“ ), 

< > 0 , 

< < ^A‘/, 

<7-/ <d\ — aj + a- + + M(1 — A-J), 

o"-/ < aj — d- + M(1 — A-J), 


□ 

Now consider the case when t) E £. Like the previous case, we need to ensure 

that after we apply the optimization strategy, t still stays the temporally closest train to the 
left of t. An analogous constraint to that of Equation (5) can be easily found by replacing t 

and t in Equation (5) with t and t respectively as follows: 

o* + d ■ - a/ - d/ 
dl + dl- 0 / - d/ 
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Note that the denominator can never be zero on the left hand side of the equation above 
because of the dehnition of t in Equation (2). To model the overlapping time ajj for all 
t) ^ £, we propose the following lemma. 

Lemma 2. For all t) E £, the overlapping time aj* between the braking phase oft on 

platform i and the aecelerating phase of t on platform j, where {i,j) G 12, ean be modelled 


by the following robust formulation 

a* — dj + e < a/ + + M(1 — A-J), (24) 

dj-a( + e<M(l-A(J), (25) 

alj > 0, (26) 

(27) 

alj < (28) 

o"-/ < dj — al + a/ + + M(1 — A(J), (29) 

< Oi — dj + M(1 — A)J), (30) 


where M is a large positive number, e is a small positive number smaller than time granularity 
eonsidered and Xj* is a binary variable whieh is one if and only if is positive. 

Proof. The lemma can be easily proved by replacing i, j, t and t with j, i, t and t respec¬ 
tively in Lemma 1. □ 

4 Modelling the constraints 

The constraints in the railway network show how the events are related. In this section we 
describe, model and justify the constraints. 

4.1 Trip time constraint 

Consider the trip of any train t G T from platform i to platform j along the track {i,j) G A*. 
The train t departs from platform i at time d\ and arrives at platform j at time a)- and the 
train can have a trip time between rb and rb. The trip time constraint can be written as 
follows: 

Vt G r V(bd) G (rb < a] -d\< rb) . (31) 
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4.2 Dwell time constraint 

When any train t E'T arrives at a platform i G it dwells there for a certain time interval 
denoted by [h*, hj so that the passengers can get off and get on the train. After the dwelling 
time is over, the train departs from the station. The difference between the departure time 
d\ and arrival time a* corresponding to the dwelling mentioned lies between h* and S^. The 
dwell time constraint can be written as follows: 

(h* <d\-a\< h*). (32) 

Every train t G T arrives at the first platform A/”* (1) in its train-path either from the depot 
or by turning around from some other line, and departs from the final platform A/”*(|A/”*|) in 
order to either return to the depot or start as a new train on another line by turning around. 
So, the train t dwells at all the platforms in A/”b This is the reason why in Equation (32) 
the platform index i is varied over all the elements of the set A/"*. 

4.3 Connection constraint: 

In many cases, a single train might not exist between the origin and desired destination 
of a passenger. To circumvent such issues, connecting trains are often used by the railway 
management at interchange stations. Let y C A/” x A/” be the set of platform pairs where 
passengers transfer between trains. If {i,j) G y, then both the platforms i and j are situated 
at the same station, and there exists a train t E T arriving at platform i and another train 
t' eT departing from platform j such that a connection time window needs to be maintained 
between train t and t' for passengers to get off from the first train and get on the latter. 
Note that order matters here. Let Cij be the set of connecting or turning-around train pairs 
for a platform pair {i,j) E y. Then the connection constraint can be written as: 

V(bi) G X V(t, t') E C,j (y“' < d^j -a\< y“'), (33) 

where y** and y() are the lower bound and upper bound of the time window to achieve the 
described connection between the associated trains. 

4.4 Turn-around constraint: 

Turn-around events provide a connection between planned terminal-to-terminal trips that 
will be executed by a single physical train. This constraint is analogous to the trip time 
constraint for these connecting movements. After a train arrives at the terminal platform 
of a train-line, it often turns around by traversing the crossing-over and starts travelling 
through another train-line. From a railway management perspective, a time window has to 
be maintained between the departure of the train from the terminal platform of the first line 
and the arrival time of it (labelled as a different train, though it is physically the same train) 
on another platform of the second line. Let cp be the set of all crossing-overs, where turn¬ 
around events occur. Consider any crossing-over {i,j) E (p, where the platforms i and j are 
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situated on different train-lines. Let Bij be the set of all train pairs involved in corresponding 
turn-around events. Let G Bij. Train t ^ T turns around at platform i by travelling 

through the crossing-over (i, j), and from platform j starts traversing a different train-line 
as train t' E T\{t}. A time window denoted by has to be maintained between the 

mentioned events, where and Kfj are the lower bonnd and upper bound respectively. We 
can write this constraint as follows: 

y{i,j) G (p V(L t') G B,j < aj - d- < Hfj). (34) 

4.5 Headway constraint: 

In any railway network, a minimum amount of time between the departnres of consecutive 
trains is always maintained. This time is called headway time. Let {i,j) G A. be the track 
between two platform i and j, and 'H^j be the set of train-pairs who move along that track 
snccessively on the order of their departures. Now, assume train t and train t' move along 
this track in same direction where G Titj. Let hf' and hf be the associated headway 
times at platform i and platform j respectively. So, the headway constraint can be written 
as: 


y{i,j) G A V(t, t') G "Hj 




tt' 


(35) 


PES2GRW2LXM2LJB2 SXZ2 ZJD2 YHR2 CSR2YSS2JYRV ^ 

□ nn o n □ nan 

□ □ □ □ CJ _n_ □ □□ □ 

YHRl CSRlYSSUYR*lS*_ 



GRW1LXML,JB1 SXZl ZJDl 


n □ 

^ LZVILHRI 


PJTIJYSILHSI 


l)ine 1 


Figure 3: Railway network considered for numerical study 


4.6 Total travel time constraint: 

To maintain the quality of service in the railway network, it is desired that for every train 
t G T, the total travel time to traverse its train-path stays within a time window [t):,,Tp], 
where and t\, are the corresponding lower and upper bound respectively. We can write 
this constraint as follows: 


yt eT (zp < aV‘(|V‘|) “ — ' VJ 


pi < 

— «V‘(1) — ^7 

Here A/'*(l) and A/'*(|A^*|) are the hrst and last platform in the train-path of t. 


(36) 
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5 Pull optimization model 

In this section we collect the objective and all the constraints discussed in the previous two 
sections, and propose our optimization problem to maximize the total duration of overlapping 
times of the SPSTPs in order to utilize regenerative braking energy produced by trains in a 
railway network. The full optimization model is as follows: 

maximize E 4‘ + E 4‘ 
subject to 

Equations (31), (32), (33), (34), (35) and (36) 
t) ^ £ Equations (5),(6)-(12) 

t ) ^ £ Equations (23),(24)-(30) 

(a->0,d^>0) 

G £ G {0, l},a“ > 0), 

where the decision variables are o*, d\, Xfj and a)*. 

As the model is a MIP with bounds, the optimization problem is AfV-haid [23, page 
242]. However, in the next section we show that for the size of the railway data considered 
in practice, the running time is quite acceptable. 

6 Numerical experiment 

In this section we apply our model to nine different problem instances to service PES2-SFM2 
of line 8 of Shanghai Metro network. The numerical study was executed on a Intel Core 
i7-46400 CPU with 8 GB RAM running Windows 8.1 Pro operating system. To model our 
problem, we have used JuMP - an open source algebraic modelling language embedded in the 
programming language Julia [24]. Within our JuMP code we have called academic version 
of Gurobi Optimizer 6.0 as the solver. In each instance we have an initial feasible timetable 
with a duration of six hours. In most of the railway networks the duration of the off-peak or 
rush hours is smaller than six hours, so a timetable that six hours is sufficient for practical 
purpose. The number of trains increases as the average headway time decreases. The results 
of the numerical study are shown in Table 1. 

The service PES2-SFM2 of line 8 of the Shanghai Metro network is shown in Figure 3. 
There are two lines in this network: Line 1 and Line 2. There are fourteen stations in 
the network denoted by all capitalized words in the hgure. Each station has two platforms 
each on different train lines, e.g., LXM is station that has two opposite platforms: LXMl 
and LXM2 on Line 1 and Line 2 respectively. The platforms are denoted by rectangles. 
The platforms indicated by PES2 and SFM2 are the turn-around points on Line 2, and the 
crossing-overs are PES2-GRW1 and LHS1-SFM2. 
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Table 1: Results of the numerical study (running time 1200 s) 


Numbei 

of 

Trains 

Numbei 

of 

Con¬ 

straints 

Binary 

Vari¬ 

ables 

Continu¬ 

ous 

Vari¬ 

ables 

CPU 

Time 

(s) 

Nodes 

Ex¬ 

plored 

Intial 

Over¬ 

lap¬ 

ping 

Time 

(s) 

Final 

Over¬ 

lap¬ 

ping 

Time 

(s) 

Initial 

Effec¬ 

tive 

Energy 

Con¬ 

sump¬ 

tion 

(kWh) 

Final 

Effec¬ 

tive 

Energy 

Con¬ 

sump¬ 

tion 

(kWh) 

Reductioi 
in Ef¬ 
fective 
Energy 
Con¬ 
sump¬ 
tion 
(%) 

100 

10437 

183 

3133 

86.64 

4534 

454.26 

2828.09 

8594.27 

7698.28 

10.42 

112 

12166 

272 

3576 

7.26 

459 

1307.23 

4354.5 

14500.05 

13740.82 

5.23 

124 

12684 

188 

3846 

0.28 

0 

1190.89 

3119.72 

8524.96 

7618.28 

10.63 

134 

13186 

128 

4081 

0.25 

0 

315.85 

1950.88 

5451.29 

4637.33 

14.93 

146 

16378 

426 

4733 

0.45 

0 

2412.92 

7314.48 

20223.63 

16998.21 

15.95 

156 

16264 

278 

4880 

15.34 

6 

2096.23 

4226.93 

14324.78 

13319.2 

7.02 

168 

18154 

390 

5346 

56.94 

1164 

1823.1 

6299.66 

20274.44 

17817.97 

12.12 

178 

19650 

472 

5723 

0.51 

0 

2874.67 

8195.41 

21685.04 

17619.35 

18.75 

190 

22947 

785 

6390 

9.6 

0 

5858.58 

12924.2, 

5 35613.96 

33574.48 

5.73 


The feasible timetables are provided to us by Thales Canada Inc.. We have taken M = 
1000 and e = 0.005. We have applied our optimization model to hnd the optimal timetable 
that maximizes the total overlapping time of the SPSTPs. We see from Table 1 that for each 
instance, our optimization model produces an optimal timetable with signihcant increase in 
the total overlapping time in comparison with the initial timetable. Such increase in the 
total overlapping time would make it possible to save signihcant amount of electrical energy 
produced by the braking trains by transferring it to the accelerating trains via the overhead 
contact lines. We can see that, in all of the cases our model has found the optimal timetables 
very quickly, the largest runtime being 86.64s. Number of nodes explored is quite small for 
every instance, and for hve out of nine instances an optimal solution is found at the root 
node. After we get the hnal timetable, we calculate the total effective energy consumption by 
all the trains involved in SPSTPs and compare it with the original timetables. The effective 
energy consumption of a train for a trip is dehned as the difference between total energy 
required to make a trip and the amount of energy that is being supplied by a braking train 
during synchronization process. So, the effective energy consumption is the energy that will 
be consumed from the electrical substations. 

The energy calculation is done using SPSIM [19] , and Cubature [25]. SPSIM is a propri¬ 
etary software owned by Thales Canada Inc. that calculates the power versus time graphs 
of all the active trains for the original and optimal timetables. Cubature is an open-source 
Julia package written by Steven G. Johnson, which is used to calculate the effective area 
under the power versus time graphs to determine 1) the total energy required by the trains 
during the trips, 2) the total transferred regenerative energy during the SPSTPs, and 3) the 
effective energy consumption that the difference of the hrst two quantities. The effective 
energy consumption of the optimal timetables in comparison with the original ones gets 
reduced for every instance, with smallest reduction being 5.23% and the largest reduction 
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being 18.74%. 


7 Conclusion 

In this paper we have presented a robust mixed integer optimization model to utilize regen¬ 
erative braking energy produced by trains in a railway network. The optimization model 
calculates a railway timetable that saves regenerative energy of braking trains by transferring 
it to suitable accelerating trains in need of energy. We have used the hypograph approach 
and interval algebra to formulate the objective function. We have modelled different con¬ 
straints to describe a feasible timetable for a railway network of arbitrary topology. We 
have applied our optimization model to different instances of a railway network for a time 
horizon spanning six hours. For each instance, our model has found an optimal timetable 
very quickly with signihcant reduction in the effective energy consumption by the trains. 
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